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ABSTRACT 

We describe the construction of a template set of spectral energy distributions (SEDs) for the esti- 
mation of photometric redshifts of luminous red galaxies (LRGs) with a Bayesian template fitting 
method. By examining the color properties of several publicly available SED sets within a redshift 
range of < z < 0.5 and comparing them to Sloan Digital Sky Survey (SDSS) Data Release 7 data, 
we show that only some of the investigated SEDs approximately match the colors of the LRG data 
throughout the redshift range, however not at the quantitative level required for precise photometric 
redshifts. This is because the SEDs of galaxies evolve with time (and redshift) and because at fixed 
redshift the LRG colors have an intrinsic spread such that they cannot be matched by one SED only. 
We generate new SEDs by superposing model SEDs of composite stellar populations with a burst 
model, allowing both components to be reddened by dust, in order to match the data in five different 
redshift bins. We select a set of SEDs which represents the LRG data in color space within five red- 
shift bins, thus defining our new SED template set for photometric redshift estimates. The results we 
get with the new template set and our Bayesian template fitting photometric redshift code (PhotoZ) 
are nearly unbiased, with a scatter of cta^ = 0.027 (including outliers), a fraction of catastrophic 
outliers (l^^phot — ^spec|/(l + ^spec) > 0.15) of 7^ = 0.12%, and a normalized median absolute rest frame 
deviation (NMAD) of (Jnmad = 1-48 x MAD = 0.017 for non-outliers. We show that templates that 
optimally describe the brightest galaxies (—24.5 < Mr < —22.7) indeed vary from z = 0.1 to z = 0.5, 
consistent with aging of the stellar population. Furthermore, we find that templates that optimally 
describe galaxies at z < 0.1 strongly differ as a function of the absolute magnitude of the galaxies, 
indicating an increase in star formation activity for less luminous galaxies. Our findings based on the 
photometry of the SDSS LRGs and our SED template fitting are supported by comparison to the 
average SDSS LRG spectra in different luminosity and redshift bins. 

Keywords: galaxies: distances and redshifts - galaxies: evolution - galaxies: fundamental parameters 
(colors) 



1. INTRODUCTION 

The primary goal of our work is to generate new spec- 
tral energy distributions (SEDs) for luminous red galax- 
ies (LRGs) to obtain precise photometric redshift mea- 
surements from template fitting methods. For the anal- 
ysis of large scale structure (LSS) based on photometric 
redshift distances it is important to precisely know the 
photometric redshift errors and biases and to make these 
as small as possible. For dark energy constraints based 
on shear tomography it is most impo rtant to know the 
errors and biases precisely ( Ma et aT|[2006| ), but on the 
other hand the rms-size of the photometric redshift er- 
rors determines how many tomographic redshift bins can 
be analyzed for a given survey. In addition, if the bias 
changes strongly as a function of the (true) redshift, then 
it cannot be treated as being constant within a given to- 
mographic interval anymore, which makes the analysis of 
such data more difficult. Another reason why small pho- 
tometric redshift uncertainties are advantageous is the 
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following. When the baryonic acoustic oscillation (BAO) 
scale is extracted from galaxy clustering measurements 
obtained w ith dista nces based on photometric redshifts 
(see Sa nchez et al.| |2011, and Bemtez et al. 2009 fo r the 
case of fhe transversal BAO scale, and Roig et aL||2009 
for the case of the radial direction), the quality of the 
result dramatically depends on the photometric redshift 
errors (see Figure 3 of|B enitez et al.||2009[ and Figure 1 
of Sanchez et al.||201ll). When galaxy properties (col- 
ors, star formation rate, dark matter halo properties) are 
studied as a function of their environment density based 
on ph oto metric reds hift distances (see, e.g., Capak et al. 
"2007' and 'Ch uter et al.||2011| for an example), redshift 
errors introduce biases in density, because galaxies are 
being preferentially scattered out of high density envi- 
ronments into the low density environments. In these 
cases it is also very important to have a low SED depen- 
dence of bias, or to at least know the bias very well as 
a function of SED and redshift. Finally, redshift errors 
imply large errors of the estimated absolute luminosities 
and stellar masses (see Nakajima et al.|2012l, Figure 13). 
If one then studies average galaxy properties as a func- 
tion of absolute luminosity or stellar mass, one has to ac- 
count for contamination from galaxies originating from 
other luminosity bins. 

In summary, it is important to know bias and redshift 
scatter, but it is also important to aim for a scatter that 
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is as small as possible, since any scatter corresponds to 
a loss of information. 

We will generate new template SEDs which optimally re- 
produce the LRG colors by fitting composite stellar pop- 
ulations (CSPs) and an extincted starburst component 
to the broadband photometric Sloan Digital Sky Survey 
(SDSS) data. We thus also obtain hints for the mech- 
anism causing the observed changes, i.e., a luminosity 
and redshift dependence of star formation history (SFH). 
Broadband photometry however is not the best method 
to study this, given that different ages and metallicities 
can produce fairly similar colors. This degeneracy can 
only be rigorously broken with spectroscopic informa- 
tion. We could use priors (i.e., exclude low metallicities) 
from the beginning, and then interpret the best fitting 
SEDs in terms of SFH. Instead, we want to make use of 
the full variety of different colors that can be generated 
even with implausible ingredients (i.e., age-metallicity 
combination, but also extinction which changes the slope 
of an SED) to increase the flexibility of SEDs that can 
be generated. We will see in Section |7.1[ however, that 
the broadband near-UV (NUV) and STUV (FUV) col- 
ors of our newly generated SEDs match the sequence 
of colors of quiescent to star-forming early-type galaxies 
(ETGs) classified by spectroscopic analyses. This means 
that generating SEDs by fitting galaxy colors with a com- 
posite population plus a burst population is not just a 
method to increase the freedom of SEDs in the u- to z- 
band (rest frame) wavelength ranges, but that it is very 
plausible that some of the LRGs do show signs of recent 
star formation. In this paper we concentrate on the es- 
timation of photometric redshifts (photo- 2:s) for spectro- 
scopically observed LRGs from the SDSS Data Release 7 



( York et al. 
consist of E'^ 



2QQQ[ [Abazajian et al.||2QQ 9). LRGs mostly 
' Gs and show rather uni form SEDs (Schnei^ 



der et al.|1983 Eisenstein et al. 2QQ3| ). Most importantly, 
their SEDs have a pronounced break at 4000 A which is 
due to an accumulation of metal lines. The 4000 A break 
together with the spectral uniformity are the reasons why 
photometric redshift estimations are very precise for red 
galaxy types. LRGs are also of great cosmological inter- 
est: They are among the most luminous galaxies exist- 
ing, are strongly correlated with each other, since they 
are preferentially located in dense regions of the universe, 
they closely map its LSS. 

This paper is organized as follows: In Section|2]we briefly 
present the photometric redshift code used in this paper. 
After that we will shortly explain the data set and inves- 
tigate the differences between the colors of the data and 
specific model SEDs (Section [3|. In Sectionjilwe present 
the method with which we generate new SETJs based on 
the LRG photometry. The photometric redshift results 
using the new template SEDs are analyzed in Section [5] 
and compared to public SDSS photo-zs in Section |6] A^ 
discussion about the properties of the new SEDs is given 
in Section [7| where we investigate the FUV and NUV 
colors in Section |7.1| and show systematic variations in 
the SEDs in Section [721 We show the results of SED 
fitting different model sets to the LRG photometry in 
Appendix [A] The colors of the created model SEDs will 
be investigated in greater detail in Appendix [B] 
Throughout this paper we assume a cosmology with 
ftm 0.3, = 0.7 and Hq 70kms"^Mpc"\ while 



magnitudes are given in the AB system. 

2. PHOTOMETRIC REDSHIFTS WITH TEMPLATE 
FITTING METHODS AND THE PhotoZ CODE 

The techniques to estimate photometric redshifts can 
be divided into two main categories: empirical and tem- 
plate fitting methods. The former use a training set of 
galaxies with known spectroscopic redshifts and photom- 
etry to derive an empirical relation between the red- 
shift and the photometric observables. The latter fit 
template SEDs to the photometry of objects, thereb}^ 
determining the most probable redshift (lAr nouts et al. 
[19991 |Be nitez|[2QQQl [Bolzonella et al ||2000 nGsabai etaT. 



2000; Bend er et al.|2001||Le Borgne Rocca-Volmerange 
2002; Babbedg e et al |2Q04| [Feldmann et al.|20061 |Bra~ 



mer et al. 2008|. The simplest approach is to perform 
solely a minimization of the difference between the 
expected colors of a SED template and observed pho- 
tometry, whereas hybrid photometric redshift methods 
combine a maximum likelihood approach together with 
an optimization routine to enhance the initial template 
set (e.g., by using a spectroscopic training set). Bayesian 
methods, on the other hand, use minimization, but 
additionally provide the opportunity to introduce red- 
shift and luminosity priors which assign probability func- 
tions to the templates, thus enabling one to lift possible 
degeneracies where models are equally probable at dif- 
ferent redshifts. The model SEDs for template matching 
methods ca n either be derived empirically (e.g., Coleman 
et al.|198Q ) , from stellar population models (e.g.,IBruzual 



& Chariot 



2003), or a combination of both. The advan- 



tage ot template fitting is that it does not only provide 
photometric redshift estimates, but also rest frame prop- 
erties of the best fit models such as the SED type and 
the absolute magnitude. Second, if good SEDs are avail- 
able, they can be used with any data set (independent of 
its photometric depth, or the filter system). Since sev- 
eral SED types and redshifts can match the colors of an 
observed object equally well, it is very important to opti- 
mize the redshift and the luminosity priors for each type, 
depending on the filter set and the data, i.e., photomet- 
ric depth and galaxy properties (which, in turn, depend 
on the galaxy selection function). 

2.1. The PhotoZ Code 

In this work we use the Baye sian template fitting pho- 
tometric redshift code PhotoZ (Bender et al.||2001 ). The 
calculation of the probabilities of the mo del- redshift com- 
binations is based on Bayes' theorem: 



P(/i|C, m) (X C{C, m|/i) • P(/i), 



(1) 



where C and m denote the colors and magnitudes of 
the photometric data, and /i are the model parameters 
z and M. C{C^m\ii) is the likelihood function which is 
proportional to exp (— x^/2). The second factor is the 
prior distribution P(/i) = P(M|T) • P(^|T), with the 
probability functions of absolute magnitude M and red- 
shift z for a template T. These priors are specified a 
priori and enable us to assign specific probability dis- 
tributions in redshift and absolute magnitude for every 
SED template. Both are parameterized as proportional 

to P{x\T) = exp (^-ln(2) (^)^), where x is the most 
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probable value ofx,x = MoTX = z. The priors were 
chosen in this functional form to be able to vary the tran- 
sitions between regions of different probabilities in such 
a way that they may be both smooth and have a box-like 
shape. In this formula the integer p controls the steep- 
ness of the probability decrease left and right of x: The 
larger the steeper the decrease (odd p values make no 
sense, of course, as the integral over the probability would 
diverge). We set the z prior to have a Gaussian shape 
(p = 2), in order to have a shallow transition from high 
to low probabilities. We will make particular use of the 
z prior later on when using SEDs specifically designed 
for different redshift ranges. For the luminosity prior we 
set X = Mr = —21, (7 = 3, and p = 6, as luminosi- 
ties Mr ^ [ —24 , —18] are very rare among LRGs (see 
also Section 



7.2 



Figure |13| ) and are thereby excluded. 
Within this interval however, the luminosity probabili- 
ties are virtually equally distributed. 
Although the LRG sample contains only objects with 
z < 0.5, the allowed photometric redshift range for all 
PhotoZ runs is [0,5.0]. In this way we can assure that 
the photometric redshift results truly are below z ~ 0.5, 
although a larger redshift range is allowed. This is im- 
portant for future work with deeper surveys where the 
redshift range is not known a priori. 
The PhotoZ code has been extensively tested and was ap- 



plied to a number of p hotometric catalogs (Drory et al. 
2QQ T]|Gabasch et al .|2004; Feulner et al. 2005; BrimiouTIe 
et al. 2008;"Gabasch et al. 2008; Brimioulle et al. 201^~Tt 
nas also been successfully applied to the Pan-STARKSl 
medium deep fields a nd is part of the P Sl Photometric 
Classification Server ( Saglia et al.||2012 ). 



3. THE DATA AND MODEL SEDS 

3.1. The LRG Sample 

The data in this paper contains objects from the spec- 
troscopic set of the SP SS Data Release 7 LRG sample 
([Eisenste in et al.||2QQl ). The selection criteria for LRGs 
in SDSS separate them from the other galaxy types by 
cuts in color space, taking their evolution with redshift 
into account. Objects passing these criteria are fiagged 
accordingly and can be obtained via the Catalog Archive 
Servei0(CAS). SDSS LRGs are color-selected in such a 
way that their SEDs have prominent 4000 A breaks. We 
remove objects from the SDSS spectroscopic LRG sam- 
ple that do not pass th e color cuts developed by Pad-^ 
manabhan et alj (2005 ), which differ niar fina lly from 



the selection criteria of Eisenstein et al. (2001). For a 
reasonable evaluation of the photometric redshift quality 
later on, we require that the spectroscopic redshift er- 
rors fulfill Azspec ^ 0.002 (the redshift resolution of the 
PhotoZ runs), and that the spectroscopic redshifts were 
determined with high confidence (> 99%). The studied 
LRG sample contains approximately 140,000 objects. 
Unless mention ed otherwise, uqriz reiev to the extinc- 
tion corrected ([Schl egel et al. | 1998| model ma gnitudes 
of SDSS in the AB system (|Oke k Gunn| 



] [1983 ). 



3.2. Colors of SED Templates versus Colors of SDSS 

LRCs 
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Figure 1. Upper panels: redshift vs. colors for the three red 
model SEDs that were previously used with the PhotoZ code. 
Lower panels: color residuals between models and data. LRGs 
are displayed in gray and their median is shown by a solid line. 
Error bars show mean data errors in five redshift bins. The dashed 
line indicates the Mannucci SO SED, whereas the colors of the two 
other red SEDs created for the HDF are shown with dotted and 
dash-dotted lines. 

In order to obtain decent photometric redshift esti- 
mates with template fitting methods, the model-set has 
to represent the data well in all available colors. In this 
section we compare the color-redshift relations for sev- 
eral red SEDs with the observed SDSS LRG colors. The 
red SEDs are taken from model sets provided in com- 
monly used template fitting photometric redshift codes. 
Additionally, we examine other SE Ds from stellar pop- 
ulatio n modeling such as t hose of [Bruzual fc Chariot 
( |2QQ3| hereafter BC03) and |Maraston et al.| ( |2QQ9[ here- 
afterM09). 

3.2.1. The PhotoZ Standard Template Set 

First, we investigate the three red template S EDs that 
have been used with the PhotoZ co de until now ([Gabasch 



eFld1|2QQ4l IFeulner et aI1[2QQ5| [Gabasch et al.||:^00^1 
Brimioulle et al. 2008)] One of the models originates 



http : //cas . sdss . org/astrodr7/en/ 



^^^^^^^^ originates 

from Mannucci et al. (2001) and describes a galaxy with 
type SO. The other two SEDs were g enerated syntheti- 
call y from single stellar populations of |Bruzual fc Char^ 
lot ( 1993|) in order to describe two slightly different sub- 
samples of red objects of the Hubble Deep Field (HDF) 
North at about z - 0.2 to 0.4. The resulting two SEDs 
are rather similar, such that one of them could, in prin- 
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ciple, be dropped in a photometric redshift analysis. In 
this case, however, photometric redshifts have shown to 
be less precise than for the case where both models have 
been used. For the success of photometric redshift esti- 
mates, it is relevant that the template SEDs cover, or at 
least span the whole color range. This is ensured by the 
Mannucci SED being too red and the other two SEDs 
being bluer regarding the 4000 A break (i.e., the g — r 
color for z < 0.35; see Figure [T]). 

Inspecting the SED colors in more detail (Figurejl]) shows 
that, indeed, none of them describes the data equally well 
at all redshifts. Those created for the HDF are too blue 
in their g — r colors for z > 0.1, and also in r — i and i — z 
throughout the whole z range. The Mannucci SO SED 
on the other hand is too red in ^ — r at z < 0.35, and in 
r — zatz>0.15. It describes the observed i — z color 
sufficiently well at z < 0.25. All three models match the 
data in 1^ — ^ within the errors, except for z ^ 0.15. 
The reason why these three SEDs give roughly correct 
redshifts for red galaxies is that the g — r color has the 
strongest redshift sensitivity. Both ellipticals predict ap- 
proximately correct values of g — r for z < 0.15, whereas 
the Mannucci SO matches the data well for z > 0.32. 
For z that lie in between those ranges, red galaxies are 
either preferentially fitted by an elliptical or the Man- 
nucci SO SED, depending on the photometric errors and 
weights applied to the i and z data. For ground based 
data that are not as well calibrated as SDSS, the i and 
z band photometric errors are often scaled up to allow 
for an imprecisely known amount of telluric absorption. 
Furthermore, if an SED template color has an offset from 
the true color which is almost constant over the whole 
redshift range (as it is the case for r—i of the two elliptical 
SEDs), then there is no contribution to a systematically 
wrong redshift estimate. 

3.2.2. Models of Bruzual & Chariot 2003 
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The CSP models of [Bruzual fc Chariot] ( |2003[ ) can be 
generated with their software galaxe^ We investigate 
SEDs from the (theoretical) BaSeL library with Padova 
1994 evol utionary t racks and an initial mass function 
(IMF) of |Chabrier| (2003) and solar metallicity. The 
models have ages from 0.01 to 13 Gyr, and exponential 
star formation rates (SFRs) with e- folding timescales r 
between 0.5 and 50 Gyr (latter introduced to account for 
a constant SFR). Furthermore, we included a model that 
has no ongoing star formation and all stars were formed 
in a single burst (henceforth referred to as r = OGyr). 
SEDs become redder in color with increasing metallici- 
ties, ages and decreasing e-folding timescales. We com- 
pare only those with the data that match their median 
best. Hence, we compute the residual between the data 
and the solar metallicity BC03 models for the u — g, 
g — r^ r — i^ and i — z colors. We weighted the residuals 
by the data errors at redshifts greater than 0.1 to avoid 
contamination by objects that do not have mostly old 
stellar populations: 
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Figure 2. Upper panels: redshift vs. colors for the three BC03 
models that match the colors of the LRG data best. Lower panels: 
color residuals between models and data. LRGs are displayed in 
gray and their median is shown by a solid line. The error bars 
show mean data errors in five redshift bins. All three models have 
solar metallicity, whereas r refers to the e-folding timescales of the 
(exponential) SFH. 

The three best fitting models have r = 1.0 Gyr and ages 
of t 13 Gyr, r = 2.0 Gyr with t = 7Gyr and t 
8 Gyr. Their colors are shown as a function of redshift in 
Figure [2] The best fitting ages and e-folding timescales 
seem realistic, and reflect that the red colors of LRGs 
originate from both shorter r and higher ages. All three 
investigated models show a reasonably good match in the 
u — g color, but are too red in ^ — r at lower redshifts, 
and too blue at higher z. Their r — i color is bluer than 
the data at very small z and deviates toward redder r — i 
in the higher z range. Finally, in i — z they are too blue 
within 0.1 < z < 0.35, and too red for small z. 

3.2.3. BPZ Models 

The ''Bayesian Photometric Re dshift'' (BPZ) cod^ 
( Benftez||2000 ; Benitez et al.||2004 ) provides a variety oi 
template SEDs. One is the niean spectrum of local ellip- 
tical galaxies by|Coleman, Wu, & Weedman (19 80[ here- 
after CWW), which was extended to the U V b}^"Be mtez] 
(2000) through linear extrapolation, and to the near-lK 
with synthetic templates of GISSE L. Extinction was ac- 
counted for following [Madau ( 1995 ). The current version 



^ ttp : //www . cida . ve/~bruzual/bc2003 or ht tp : //www . iap7| 
|fr/~cha rlot/bc2003^ ' ' 



^vl.98b: http : / / acs . pha. jhu . edu/~txitxo/^^ 

'^vl.99.3: http: //www. its . caltech. edu/'-coe/BPZ/ 
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Figure 3. Upper panels: redshift vs. colors for the red CWW 
models: the original SED (dashed), the repaired SED of BPZ (dot- 
ted), and the LePhare CWW model (dash-dotted line). Lower pan- 
els: redshift vs. color residuals between models and data. LRGs 
are displayed in gray and their median is shown by a solid line. 
The error bars show the mean data errors in five redshift bins. 

of BPZ (B enitez et al.|2QQ4 ) includes a repaired version of 
this elliptical SED, which was calibrated with an array 
of ground based data. The colors of both the original 
and the repaired red CWW SED, are shown in Figure 3j 
Both match the data well in u — g and r — with the 
exception of the latter for higher z. There are deviations 
from the data in ^ — r over the whole considered z range, 
except for z > 0.45 for the original CWW SED, and for 
medium redshifts for the repaired version. Furthermore, 
they are able to match the data median in i — z only for 
medium and high redshifts. 



3.2.4. LePhare Models 
Another Bayesian photo- z sof tware 



is LePhar^ 

( [Arnouts et al.||l999[ [Ilbe rt et al. 2006), which uses a 
wide range ot template SEDs. We adopt the nomencla- 
ture of the models used within the LePhare code in the 
following. Again, we investigate only red SEDs that pro- 
duce the smallest residuals from the data according to 
Eq. [2] The SED delivering the lowest residual originates 
from CWW models. It was generated by a linear interpo- 
lation of the four CWW SEDs and a star-forming galaxy 
that was constructed with the population synthesis code 
of Bruzual & Chariot. The resulting colors of the inter- 
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Figure 4. Upper panels: redshift vs. colors of the three models 
of LePhare that match the data best: the COSMOS (dashed), the 
GRASIL (dotted), and the GISSEL SED (dash-dotted hne). Lower 
panels: redshift vs. color residuals between models and data. 
LRGs are displayed in gray and their median is shown by a solid 
line. The error bars show the mean data errors in five redshift bins. 
A description of the model SEDs is given in the text and references 
therein. 

polated red CWW SED are shown in Figurejs] together 
with the colors of the original red elliptical CWW SED 
and its repaired version of BPZ. Its colors are very similar 
to the original model and the repaired BPZ version, which 
were investigated in the last section (Figure [3|. The fact 
that these observationally created SEDs can reproduce 
the LRG colors well, and are therefore in principle well 
suited for photo-z estimations, was also pointed out by 
Benitez (2000). 

Figure [4] shows the colors of three other SEDs that yield 
the second, third and fourth smallest residuals from the 
data. The one termed COSMOS is from the LePhare set 
called COSMOS-SED. Models from this set were used in 



^ http : //www . cf ht . hawaii . edu/~arnouts/lephare . html 



Ilbert et al. (2009) for photometric redshift estimates in 
the 2 deg^ COSMOS field. The model tagged grasil is 
a 13 Gyr old elliptical galaxy t hat was created with the 
GRASIL code (|Silva et al.||1998|). The SED describing the 
data fourth best comes from the LePhare 42-GISSEL 
set. It has an e-folding timescale of 1 Gyr, a metallicity 
of Z = 0.02 (Zq) and an age of 7 Gyr. From Figure [4] 
one can infer that the g — r colors of those three models 
are redder than the data for low z. In r — z, the three 
model SEDs are too blue for lower redshifts. Their i — z 
is redder than the data at z < 0.07 and z > 0.45, and 
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Figure 5. Upper panels: redshift vs. colors for three M09 models, 
where the LRG with an age of 6 Gyr (dotted line) matches the data 
best regarding the residual to the median. Lower panels: redshift 
vs. color residuals between models and data. LRGs are displayed 
in gray and their median is shown by a solid line. Error bars show 
the mean data errors in five redshift bins. The other lines display 
the M09 models that have an age of 4 Gyr (dashed), and 12 Gyr 
(dash-dotted line). The solid circles show a 12 Gyr (at z = 0) old 
M09 SED model evolving with time. The following dots represent 
the same star formation model, but with ages of 11, 10, 9, 8, and 
7 Gyr. A description of the model SEDs is given in the text and 
references therein. 

bluer at 0.07 < ^ < 0.3 

3.2.5. Models of Maraston et al (2009) 

Maraston et al. (|2009|) already addressed the problem 
of the color mismatch between model SEDs for LRGs and 
the actual data. They studied several SEDs and com- 
pared them to the data median, focusing on ^ — r and 
r — i, and discovered that neither the synthetic nor the 
empirical SEDs could accurately represent the data at 
all redshifts equally well. To produce new models that 
resemble the color median of the data more precisely, 
Maraston et al. added a mass fraction of 3% blue metal- 
poor stars to create empirically corrected CSP models. 
It enabled them to match the data colors considered in 
their publication {g — r and r — i) very well. We present 
the colors of a subset of the M09 model SEDqf] in the 
u — 9^ 9 — ^ — h and i — z colors in Figure [5] The mod- 
els that are displayed by lines have not been evolved with 

^ i h ttp: //www. icg. port . ac. uk/~maraston/Clau(iia' s_ | 
|Stellar_Population_Model . htmTl 



look-back time. The SED with an age of 6 Gyr yields the 
smallest residual (according to Equation [2| to the data. 
It increases faster when turning to younger populations 
than it does for older ages. In order to show the behav- 
ior for ages different from 6 Gyr, we display one older 
(12 Gyr) SED in Figure [sj and one younger model with 
an age of 4 Gyr. 

In contrast to the SEDs examined above, the M09 red- 
shift-color relation has the same shape as that of the 
data median. The u — g color deviation of the M09 mod- 
els from the data lies well below the mean error, and 
does not differ notably for the models with 6 Gyr and 
12 Gyr. At z < 0.2 the u - g color of the 4 Gyr SED 
resembles that of the other two SEDs but has lower val- 
ues for higher redshifts. The g — r color of the 6 Gyr old 
LRG is close to the data at low z {< 0.25) but deviates 
at z > 0.35. Its r — z and i — z colors are too blue in 
comparison with the data at z < 0.4 and z < 0.3 respec- 
tively. The 12 Gyr old SED matches the data in r — z 
well for z < 0.4 and in i — z for all whereas its ^ — r 
color is too red in the whole redshift range. The LRG 
with an age of 4 Gyr is bluer with respect to the data 
in r — z at all redshifts and i — z at z < 0.3, as well as 
in ^ — r at intermediate redshifts (0.2 ^ z ^ 0.35). Un- 
surprisingly, the older SEDs represent the data at lower 
redshifts better, while the younger population exhibits 
colors matching better at higher z. An exception is — r 
where the oldest model shows an offset from the data, 
whereas the 6 Gyr old model fits well up to z ^ 0.35. 
The SED with an age of 4 Gyr matches the data reason- 
ably well in ^ — r throughout the whole redshift range. 
Figure |5] shows furthermore the M09 model with an age of 
12 Gyr at z = evolving with look-back time. The dots 
at higher z represent ages of 11, 10, 9, 8, and 7 Gyr. The 
colors of the redshift-evolved M09 model are in very good 
agreement with the data compared to the previously in- 
vestigated model SEDs. However, the model which is 
12 Gyr old at z = and evolves with redshift exhibits an 
overall offset from the data median in g — and has a 
bluer r — i color for 0.15 < z < 0.35. Furthermore, it 
deviates in — ^ at higher redshifts {z > 0.3), but still 
lies well within the errors. 

3.2.6. Summary 

Altogether, we conclude that we could not find models 
that are able to represent the SDSS LRG data within 
their errors at all redshifts. The reason for this is that 
the SEDs differ from the data median as a function of 
z. The M09 models (not taking evolution with look-back 
time into account; shown with lines in Figure]5| describe 
the SDSS LRGs best compared to the templates exam- 
ined before, but they also suffer from notable deviations 
within some redshift intervals. The lower redshift range 
(z < 0.1) and the i — z colors were not considered in 



Maraston et al | ( |2009 ). 

As expected, the colors of the redshift evolved SED with 
an age of 12 Gyr at z = (solid circles in Figure [5|) 
match the data better as a function of redshift, although 
not perfectly. Furthermore, if only the age of a SED 
is varied, one just obtains a one-dimensional sequence 
in color-color space. Instead, the colors of true SEDs 
are spread in two dimensions, which is a clear hint that 
not only the age of the LRGs, but also the SFH, and 
maybe metallicity, vary. A wider spread can thus only 
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be produced if different SFHs and metallicities are con- 
sidered, which are not pubhcly available for M09 models. 
We therefore conclude that in order to get SEDs which 
represent the data at all redshifts and in all colors, one 
has to consider different SED template sets for different 
redshift bins. 

4. NEW SED TEMPLATES 

In the past, optimal template SEDs have been ob- 
tained by "repairing" individual SEDs empirically using 
the mism atch of observed and model colors (e.g., Csabai 
|et al. 2003). In addition, or alternatively, a larger variety 
of template SEDs has been obtained by interp olating be 
tween a sm aller number of template SEDs (e.g.. 



Brammer 



et al.|[2QQ8 ). However, there are some caveats. If models 



that one starts with are a bad match to the data it is 
unlikely that a combination of them will be any better. 
Moreover, if one introduces SEDs that do not match the 
data, one deteriorates the photometric redshift quality. 
Furthermore, the question of how many SEDs are needed 
to match the data remains undecided. 

4.1. Generating Best Fitting SEDs for Individual 
Objects in the LRG Catalog by SED Fitting 

We start with SEDs that match the photometry of in- 
dividual galaxies as well as possible, and subsequently 
select a small subset of these SEDs with the goal that 
they together describe the whole data set. To derive 
SEDs for individual galaxies, we make use of the largest 
appropriate freedom to describe a SED with five band 
data. I.e., we describe the SED of an LRG by a super- 
position of a CSP and a burst model. Furthermore, we 
allow for extinction of the CSP and the burst component. 
The burst and its extinction are physically motivated by 
the fact that LRGs can contain some young stars. The 
extinction of the CSP is mainly there to increase the 
degree of freedom and allows to change the continuum 
slope (i.e., to introduce a "variance" of CSP models) of 
the main stellar population. 

We generate SEDs for LRGs with the SED fitting routine 
SEDf it (Dro ry et al.||2QQ4| ). It chooses the best fitting 
CSP+burst combination through minimization, con- 
currently allowing for dust red dening of both compon ents 



following the extinction law of Calzetti et al. (12000) 



In Section [3?2] we pointed out that for each model there is 
always a discrepancy between the model and data colors 
at some and we are not able to pick a set of models 
which describe the data equally well at all z. Thus, we 
create templates specifically designed to represent dif- 
ferent redshift regions. We do so by splitting the spec- 
troscopic SDSS LRG sample into five subcatalogs, all 
of them containing objects within redshift intervals of 
width 0.04, centered on z = 0.02, 0.1, 0.2, 0.3 and 0.4. 
These catalogs contain 9660, 38,583, 9292, 15,093, and 
11,667 objects. They are then divided into two equally 
sized subcatalogs; one serving for template generation 
(derivation half) and the other one serving as reference 
when estimating photometric redshift accuracies (valida- 
tion half). Subsequently, we run SEDf it with a model set 
consisting of 864 BC03 m odels from the BaSeL library 
dBruzual fc Charlot|20 Q3[) with Padova 1 994 evolutionary 
tracks and an IMF from |Chabrier ( 2003 ) . They have ages 
between 0.01 and 13 Gyr, e- folding timescales of r = 0.0, 
and r = 0.5 to 50 Gyr, and four different metallicities. 



Z = 0.004, 0.008, 0.02 (Zq) and 0.05. In order to fit the 
photometry, these models can be superposed by a frac- 
tion of the burst model, with an essentially constant SFR 
(i.e., r = 20 Gyr), solar metallicity, and an age of 50 Myr. 
The burst mass-fraction is constrained to [0.0, 0.01] with 
steps of 0.002. Both the CSP and the burst component 
can be extincted separately. This fitting procedure re- 
sults in a SED that reproduces the data as well as pos- 
sible. We end up with one "best fitting SED" for each 
object in the template generation catalog. 
The relative performance of SED fitting with models of 
BC03, M09, and Maraston & Strom backl|201l[ hereafter 
Mil) is demonstrated in Appendix [a{ wnere we also jus- 
tify why we prefer BC03 models overM09/ll for our SED 
construction procedure. From Figure [T9| in Appendix [A] 
we conclude that BC03 models allow a better coverage 
of the colors of the LRGs than M09/11 models do, and 
fit the data better (Figure [20]). 

4.2. Selection of Best Fitting SEDs for the New 
Template Set 

We want to select a set of models from all the SEDs we 
produced for the template generation subcatalogs such 
that this set represents the five redshift subcatalogs in 
color space. We explain our procedure for the z ~ 
redshift interval using Figure |6] 

First we plot the distributions of the measured colors 
of the galaxies in the u — g versus r — i plane. Objects 
which have a u band error smaller (larger) than the me- 
dian are displayed with light (dark) blue points. Second, 
we plot the colors of the best fitting SED for each ob- 
ject with red points. If the best fitting SEDs are a good 
match to the galaxy data, then the distribution of them 
in color space should be similar to the observed distribu- 
tion, once the spreading by photometric errors is taken 
into account. The mean photometric errors for the two 
subsamples (with high and low u band errors) are shown 
in dark and light blue in the lower right corner of the 
panels in Figures |6j and 21- 25 It shows that galaxies 
with larger photometric errors in u in general also have 
larger photometric errors in r — i and i — z, which im- 
plies that these objects are fainter in all observed bands, 
rather than just in the u band. We conclude from the 
shown sizes of the mean photometric errors on one hand 
that the spread in colors is of physical origin, and not 
due to photometric errors. This also holds for higher 
redshifts where the error bars are of significant size. In 
order to confirm that, we simulated galaxy colors by as- 
suming that the spread in color space is a result of the 
photometric errors only. We simulate the data by as- 
suming that each object in fact has a color value of the 
average color. Then we assign a random deviation from 
this point, based on a Gaussian probability distribution 
with the measured error as standard deviation. Since the 
spread in the simulated data is smaller than that of the 
data, we confirm that the spread within the colors is of 
physical origin. On the other hand, we infer that the 
color distribution should hardly be broadened by photo- 
metric errors, and that the best fit SED color reproduce 
the variety of 'true' galaxy colors well. Only for the red- 
dest colors {u — g > 1.5 and r — i > 0.45) we find mostly 
objects with larger than median photometric errors, such 
that some of these data points will be scattered out from 
the true distribution by large photometric errors. 
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Figure 6. Color vs. color plots for SDSS data (blue) and their 
individual best fitting SEDf it-SEDs (red) for 0.0 < z < 0.04. Ob- 
jects with u band errors lower than the median are indicated in 
light blue, whereas those with greater errors are dark blue. The 
grid (upper panel) and the boundary (lower panel) within which 
the models are selected are shown in black. The dark green dots 
are the preselected models, whereas the light green points repre- 
sent the models that are left over after the removal of redundant 
SEDs. The yellow dots are models that shall account for objects 
outside the selected boundaries in ugri and griz. 

We now select model SEDs for the "main galaxy" pop- 
ulation. Therefore, we define an area in the ii^ri-plane 
which contains the majority of objects (black box), and 
then set a grid within this area (black grid). The cho- 
sen area is selected by eye and was confirmed to include 
at least 70% of the objects from that redshift bin. We 
want to select one object per cell which shall represent 
the other objects within that cell. We carry out this pro- 
cedure only for cells that contain at least 20 galaxies. 
From each of those cells we pick the SEDs best fitting 
the five objects from that cell which have the lowest u 
band error. 

As mentioned above, our objective is to produce pho- 
tometric redshift s with biases as small as possible, and 



we therefore try to single out SEDs that produce small 
(Az) on the catalogs of the regarded z range. Hence, 
we perform PhotoZ runs on the derivation half and on 
the validation half of the z ~ 0.02 subcatalog, the lat- 
ter to confirm that we did not unintentionally divide the 
catalog into halves with different object properties. We 
perform the photometric redshift estimation with each of 
the preselected models separately, thus fitting only one 
model per run. Afterward, we pick the SED that renders 
the smallest bias on both subcatalog halves out of the 
five originating from the same cell. The selected model 
shall represent the objects within its neighborhood in 
color space in the following PhotoZ runs. Repeating the 
same procedure for all cells and redshift bins, we select 
roughly 20 to 35 model SEDs per z bin. 
In order to save computation time for a PhotoZ run, we 
want to keep the number of templates to a minimum. 
Therefore, we calculate the root mean squared fiux dif- 
ference of each SED combination at every defined wave- 
length. The most similar SEDs are then removed from 
the final set, leaving typically about 10 models per red- 
shift bin with a total number of 49 templates. 
The redshift of the models is plotted versus their color 
in Figure [7| Color-color relations of the new models in 
comparison to the LRGs are displayed in Figure |6| and 
Appendix |Bj Figures 21-25] The templates represent the 
LRG colors and their spread for most values of z. The 
exceptions at z ~ 0.4 are discussed in Appendix [Aj Fig- 
ure [25] We test the performance of the new models with 
respect to photometric redshifts in the next section. 
As an alternative we also selected template SEDs with 
a kd-tree. We generated one hyperplane for each of the 
four colors u — 9 — r — and i — z^ thus a four- 
dimensional tree with 16 template SEDs for each red- 
shift bin. One model was selected from each branch of 
the tree by the same criteria as already described. As 
this technique yielded SEDs not too much different, we 
kept our hand-selected set. 



5. 



PHOTOMETRIC REDSHIFT PRECISION WITH THE 
NOVEL TEMPLATE SEDS 



Before analyzing the photometric redshift results we 
get with the new template SEDs, we introduce some 
quality parameters. To define badly estimated photo- 
metric redshifts we use the threshold l^^nho t 



^spec) > 0.15 introduced by Ilbert et al. (2006). Ob- 



^sp ec I (1 ~l~ 



jects that fulfill this inequality are tagged catastrophic 
outliers", and their number fraction will be called 77. 
aAz shall denote the root mean square of the photo- 



metric redshift error Az = 



^phot 



A parame- 



ter that quantifies the distribution perpendicular to the 
^spec = ^phot line without the inclusion of catastrophic 

outliers is (Jnmad = 1.48 x median (^j^ 



\Azl 



non— outliers 

dllbert et al.|[2Q06l ). With cnmad one can make conclu- 
sions about the redshift uncertainty excluding the tails 
of the distribution. 

After these definitions we turn to the integration of the 
newly created model SEDs into PhotoZ. In Section [2] we 
mentioned the redshift priors we can impose on each tem- 
plate. The total probability for a model-z combination is 
proportional to the redshift probability P^. As already 
mentioned in Section [2j we use a Gaussian probability 

distribution, Pz{z) ^ exp ^— ln2 In this way, z 



Photo-zs of LRGs from the SDSS DR7 



9 




0.0 
0.5 r 



2.0 



0.5 
0.2 

S 0.0 
9^-0.1 
-0.2 



: 1 1 1 1 : 

— » 0) $ 1 { ~ 


^ 1 1 1 1 

-a i 1 J 1 - 


r 1 1 1 1 


1 1 1 1 ^ 

h \ ' i ' It ' I i 


i! i ■ 1 M M : 

r • • "i 

\» ' 1 1 1 — 


r — 1 1 1 1 — i 



0.0 



0.2 
0.1 



0.0 

-0.1 

-0.2 



0.0 0.1 0.2 0.3 0.4 0.0 0.1 0.2 0.3 0.4 
redshift redshift 

Figure 7. Upper panels: colors of the novel template SEDs as 
a function of redshift. Lower panels: redshift vs. color residuals 
between models and data. LRGs are represented by the density 
plot, whereas the solid line is their median. The error bars show 
the mean data errors in five redshift bins. The colors and the 
color residuals of the new template SEDs are represented by points. 
They are plotted at the mean values of the redshift bins from which 
they were created. 

can be set to the according redshift interval's center for 
each model, thus z = 0.02,0.1,0.2,0.3,0.4. The redshift 
prior probability of the model is equal to ^ at the corre- 
sponding z ±0.2. 

The photometric redshift results are shown in the upper 
panel of Figure [Sj With the new red templates, the de- 
scribed priors, and a resolution of 0.002 we indeed get 
very accurate redshifts. The catastrophic outlier rate is 
very small at = 0.12 %. The dispersions is ctaz = 0.027 
and ctnmad = 0.017 respectively. The mean absolute de- 
viation of Zphot from Zspec is (|Az|/(l + 2;spec)) = 0.015. 
We show the bias as a function of redshift in the middle 
panel of Figure [Sj It exceeds 0.01 only for the highest 
redshifts {z > 0.47). The reason for the small positive 
bias at z < 0.02 is that there are no negative photometric 
redshifts that can cancel out overestimation. The total 
mean bias reads -0.0004. At z^^^c ^ 0.36, the 4000 A 
break position in wavelength cannot be determined to 
high precision. This is because at z ~ 0.36 the 4000 A 
break is observed at A ~ 5500 A and thus lies in the gap 
between the SDSS g and r filter. We therefore expect 
that the error distribution widens at this redshift. The 
reason for the error becoming bimodal is likely because 
at this redshift we do not have "perfect templates" and 
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Figure 8. Photometric redshift results of a PhotoZ run with the 
newly created templates and 140,331 SDSS LRGs. We removed 20 
objects that had undefined photometric redshifts. ctaz^ ctnmad^ 
and T] are explained in Section [s] 

thus the photometric redshifts are either over- or under- 
estimated. Luckily the errors nearly average out for every 
redshift bin. 

The bottom panel of Figure |8] shows the dependence of 
the scatter aAz and ctnmad on z. At low z (< 0.01) 
it is up to 0.04, which is likely due to contamination 
of the LRG sample with non LRGs. This is supported 
by the fact that the values in this region are signifi- 
cantly higher than for z > 0.01. Apart from very small 
redshifts, the scatter stays more or less the same for all 
redshifts. 

6. COMPARISON TO SDSS DATABASE RESULTS 

We now compare our photometric redshift results with 
those in the SDSS database. There are three kinds of 
photometric redshifts available in the database: One was 
obtained through template fitting (i.e., a hybrid tech- 
niq ue of template fittin g and a reparation of templates, 



Csabai et al. 2003|, and the other two are results 



of an a rtificial neural network (ANN) approach (Oyaizu 



et al.| 2QQ8). The results of the template fitting are shown 
in Figure [9) On average, the a^z and (Jnmad of these 
redshifts are smaller than for our method, but the bias 
is more unevenly spread in redshift space. Also for this 
method the median bias is smaller than 0.01 only in the 
redshift range of 0.07 < z < 0.32. For other redshifts the 
bias exceeds that of our redshifts considerably. 
The SDSS database values for photometric redshifts ob- 
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Figure 9. Photometric redshift results of SDSS template fitting. 
We removed those objects with undefined photometric redshift s. 
V: ctnmaDj and (Jaz cire explained in Section [s] 
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{CC2), re- 



tained with the ANN are shown in Figures 10 and [TT| 
where we display the results of two techniques, Dl and 
CC2. The ANN(Dl) method uses ugriz magnitudes and 
the concentration indices of all five bands for the determi- 
nation of photometric redshifts. A concentration index 
of a passband i is defined as the value of the ratio be- 
tween the radii encircling 50% and 90% of the Petrosian 
fiux. The CC2 approach on the other hand is based on 
the u — g — f'i r — i^ and i — z colors and the concentra- 
tion indices in the r and i bands. Photometric redshift 
results are shown in Figures [To] {Dl ) and 
spectively. 

The ANN(Dl) performs excellently in terms of scatter 
and bias from z = 0.01 to z = 0.3, but it yields a 
considerably larger bias than our redshifts for z > 0.3. 
Photometric redshifts from the CC2 method are signifi- 
cantly overestimated for z < 0.08, where the bias exceeds 
0.01. Additionally, the ctnmad is particularly high at 
small spectroscopic redshifts {z < 0.06), which is proba- 
bly due to the overestimation in this range. At z ^ 0.3, 
Az surpasses 0.01, whereas at intermediate redshifts, the 
photometric redshifts from the ANN(CC2) are slightly 
overestimated. 

We conclude that PhotoZ with the new set of LRG 
templates delivers results competitive with empirical ap- 
proaches. 

7. PROPERTIES OF THE NOVEL TEMPLATE SET 
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Figure 10. Photometric redshift results of SDSS ANN(Dl). We 
removed those objects with undefined photometric redshifts. 77, 
c^NMAD, and (j/\z are explained in Section [s] 

In this section we want to discuss the properties of 
our new template SEDs in the UV wavelength range and 
analyze how well they represent objects with UV data. 
Being able to match the UV colors of observed objects 
justifies that we discuss the novel SEDs from a physical 
point of view and wit h re gard to the underlying galaxy 
population in Section [7^ 



7.1. UV colors of the Novel Templates and MOO Models 

We show color-magnitude and the color-color relations 
of our final template SED set down to UV wavelengths 
within 0.05 < z^^^c < 0.1 2 and compare our findings 



to those of Ree et al. ( |2Q12 , hereafter R12). We want 
to evaluate if the novel models are able to represent 
the Galaxy Evolution Explorer (GALEX) data, since we 
used no information about the UV in the creation of our 
SEDs, as even the GALEX NUV filter central rest frame 
wavelength is Anuv, central = 2544 A at a redshift of 0.12, 
whereas the SDSS u band filter has a much larger central 

wavelength of A^, central = 3546 A. 

R12 selected ETGs in the SDSS data on the basis of 
morphology, requiring f racDeV_r ^ 1, and removed con- 
tamination by late-types through visual inspection. The 
sample galaxies were then classified into quiescent, star- 
forming, composite, and active galactic nucleus (AGN) 
categories on the basis of their emission line characteris- 
tics. Quiescent ETGs were selected in such a way that 
they do not exhibit emission lines. 
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Figure 11. Photometric redshift results of SDSS ANN(CC2). ry, 
(Jnmad? cind a/\z are explained in Section [s] 

Figure [12] shows the data medians and their spread of 
the four ETG categories of the galaxy sample of R12 in 
color-magnitude {u — r versus Mn and FUV — r versus 
Mr) and color-color diagrams {u — r versus u — g and 
FUV — r versus FUV — NUV) in analogy to Figure 1 in 
R12. Additionally, we plot the colors and magnitudes of 
our final templates after fitting them to the spectroscopic 
redshift of the LRGs within 0.05 < z < 0.12, as well as 
their intrinsic colors within that z range. The median 
values of u — g and u — r coincide with those of the qui- 
escent and AGN ETGs of R12 (lower left panel), and we 
are able reproduce the observed UV-colors. med(Mi^) 
on the other hand is about 0.4 mag lower in this work 
when compared to the quiescent and AGN ETGs of R12 
(upper panels), visualizing that the objects in our sam- 
ple are definitively LRGs. Furthermore, the median UV 
excess (FUV — r) of SEDs of this work is significantly 
smaller than that of R12 (~ 1 mag compared to the qui- 
escent galaxies; upper right panel). The lower right panel 
displays the UV-excess versus FUV — NUV, as well as 
the UV red sequence found by R12 (solid line) which is 
derived empirically through a linear fit to the observed 
colors with width 3 x MAD (dotted lines). A fraction of 
our models however lie above the UV red sequence. Also 
the FUV — NUV color is higher than in R12, due to the 
fact that our models have higher fluxes in the NUV band 
and lower fluxes in FUV. According to [Schawinski et aL 
|2007), recently star forming (RSF) galaxies fulfill the 
semi-empirical criterion NUV — r < 5.4, which is plotted 



in Figure [12] (dashed line). A considerable part of our 
models is tlius classified as RSF. 

The lower panels of Figure 12 show M09 models on 
top of the SEDs of this work. In comparison to R12, 
M09 models neither cover the region of higher UV-excess 
(FUV - r < 6), nor those where FUV - NUV < 1.7, and 
n early all o f them exhibit no RSF following the criterion 
of |Schawmski et aL, (|2007j) . M09 models account for pas- 
sive evolution only, which is why they do not populate 
these RSF regions. However, the M09 SEDs match the 
u — r versus u — g colors of quiescent, AGN and com- 
posite ETGs from R12. Nevertheless, at low redshifts 
the LRG sample consists not only of passively evolving 
galaxies, but shows signs of recent star formation. Mil, 
which have a large variety of SFHs, ages, and metallici- 
ties, should be able to cover these regions. 

7.2. Differences in the SEDs within z — Mr Bins 

By construction our new template SEDs do not neces- 
sarily consist only of old CSPs (which are necessary to 
match the 4000 A break of the LRGs), but can have im- 
prints of young stellar populations (i.e., more UV-excess 
relative to old populations) and dust extinction, chang- 
ing the slope of the CSP and the burst. Instead of ex- 
amining what kind of combination the LRGs are fitted 
by, we present how different the SEDs of galaxies are 
when matched by our templates. Therefore, we deter- 
mine the best fitting SED for every LRG by fitting our 
new SED set at the spectroscopic redshifts of the objects 
and calculating the absolute magnitude in the R band 
filter thereafter. We then bin objects within the z versus 
Mr plane (Figure 13) and calculate the mean SED of 
that bin, by adding up all best fitting SED templates of 
objects from that bin with weights given by their flux in 
R. Doing this for different redshift and luminosity bins, 
we obtain the luminosity weighted light distribution of 
all objects within a given redshift and luminosity range, 
under the assumption that they are well described by one 
of the SEDs from our template set. The results for lu- 
minosity bins with z < 0.1 are shown in Figure [14] The 
SEDs which match the SDSS data best show a strong 
dependence on the brightness of the object: The fainter 
galaxies are, the more flux they have at A < 3500 A rela- 
tive to the brightest Mr ^ —24 bin, and the steeper the 
slope of the SED gets for A > 5000 A on average. Both 
observations can be explained by star formation in these 
galaxies, i.e., that low mass LRGs still build stars at low 
z in contrast to the more massive LRGs. Hence, the di- 
versification in the u — g color is not reached by simply 
selecting smoothly varying SFHs, but rather by adding 
to these single starburst events of differe nt amounts. We 
will confirm this result in Section [7.2.11 
Figure [15] shows the superposition of the models for 
z > 0.1 and —24.5 < Mr < —22.7 in six consecutive red- 
shift bins. At blue wavelengths (A < 3000 A) the SEDs 
from objects at higher redshifts exhibit more flux than 
those at lower redshifts, which indicates increasing star 
formation turning on to higher z. These effects are in ac- 
cordance to what one would expect for older stellar pop- 
ulations at low z and younger objects at higher z. Older 
red populations have (nearly) completed star formation 
and hence show low fluxes in bluer bands, whereas their 
fluxes in the red (A > 6000 A) exhibit a smaller contin- 
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Figure 12. UV colors of the new template SEDs from a fit to the spectroscopic redshift of the LRGs within 0.05 < 2;spec ^ 0.12 
(gray/black). Yellow filled circles with error bars indicate the median colors and magnitudes with standard deviations. Medians and 
spreads of R12 are displayed in red (quiescent), green (AGN), cyan (composite) and blue (star-forming ETGs). The orange sequences of 
dots joining a line segment in the lower two panels indicate the new models from SED fitting, whereas those in purple illustrate the colors 
of M09 models, both within the same redshift range. The solid line represents the UV red sequence of R12, the dotted lines are 3 x MAD. 
The criterion of RSF galaxies (NUV — r < 5.4) found by Schawinski et al. (2007) is displayed by a dashed line. Plotting ranges correspond 
to those of R12 (Figure 1). 



uum slope compared to younger galaxies. The SED from 
the first redshift bin (0.1 < ^spec < 0.17) is slightly bluer 
than the following (0.17 < ^^gpec < 0.23), because there 
are still star forming objects within that bin. This can be 
seen in Figure [Tsj there are still fainter galaxies in this 
bin, wherefore not all of them a re pas sively evolving. Ad- 
ditionally, Figure [16] in Section [7. 2 . 1 1 confirms that there 
are still star forming galaxies in this 2:/Mi^-range. 
In order to support the latter discussion quantitatively, 
we performed a to the variety of BC03 models 

that are the basis of our template generation in Sec- 
We keep the metallicity fixed at Z = as 
= 0.0 Gyr in order to avoid degeneracies be- 
r, a nd ag e. (We discussed these degenera- 



tion 4.1 



well as r 
tween Z, 

cies in Section |3.2.2[ ) The best fitting ages for the red- 
shift-luminosity bins at z < 0.1 are 9, 9, 7, 4.5, 4.5, and 
3.5 Gyr from brighter to fainter bins. Keeping only the 
metallicity fixed at Z = Z©, we observe that the best 
fitting models nearly all have r = 1.0 Gyr, which sup- 
ports the claim that they all exhibit ongoing star forma- 
tion. The first (most luminous) bin has a best fitting r of 
0.0 Gyr, which is due to the very red SED we calculated 



for it. Anyway, since this bin contains only one object we 
cannot assume the SED to be representative, although it 
would be in good agreement with the sequence the su- 
perpositioned SEDs form, and it w ould n ot contradict 
the observation we make in Section 17.2.11 where we de- 
duce from the relations between emission lines that no 
star- forming galaxies lie in this particular z — Mr range. 
We also performed the fitting of BC03 models to the 
superpositioned SEDs for 2: > 0.1 and —24.5 < Mr < 
—22.7 keeping Z at Zq and r = 0.0 Gyr fixed. The best 
fitting ages are then: 12, 13, 11, 9, 9, and 9 Gyr with 
increasing z. As was mentioned above, the first redshift 
bin (0.1 < Zgpec < 0.17) still contains objects with ongo- 
ing star formation, which is why the best fitting age is 
lower than in the second bin (0.17 < ^^gpec < 0.23). Nev- 
ertheless, the ages from the other five redshift slices are 
monotonically decreasing. By keeping only the metallic- 
ity fixed, the r values increase from lower to higher red- 
shift, which also yields to bluer spectra. Furthermore, 
the differences in the light travel time between the cen- 
ters of the first and the last redshift bins is ~ 3.1 Gyr, 
which is in agreement with the results, and also 
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supports what we expect when we assume that the un- 
derlying population of our sample at z > 0.1 predomi- 
nantly consists of passively evolving galaxies, and that 
the variations in the SEDs are a result of aging. 
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Figure 13. Density plot of spectroscopic z vs. absolute R band 
magnitude. is obtained from the template SED which matches 
a galaxy best at its spectroscopic redshift. The solid lines show bins 
within z or Mji intervals. The luminosity weighted means of the 
template SEDs that match each object in a given z — Mr bin best 
are shown in Figures |l4| and [Ts] 
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Figure 14. Flux weighted superposition of model SEDs (normal- 
ized to fg) for objects with 0.0 < Zspec < 0.1 within consecutive 
luminosity bins. The solid gray lines are the normalized SDSS fil- 
ter curves, whereas the dashed gray lines represent the filter curves 
sX z = 0.1. The bins contain 1, 6943, 23,501, 1757, 123, and 33 
objects (from higher to lower luminosities). 



7.2.1. Comparison to SDSS Star Forming LRGs 

We want to compare our findings of the previous chap- 
ter to the observed spectra of SDSS LRGs. There- 
fore, we download the SDSS spectral subclassifications 
of our LRG sample that were produced by the SpecBS 
pipeline, where galaxies are further categorized by their 
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Figure 15. Flux weighted superposition of model SEDs (normal- 
ized to fg) for objects with —24.5 < Mji < —22.7 within consec- 
utive redshift bins. The solid gray lines are the normalized SDSS 
filter curves at z = 0.1, whereas the dashed gray lines represent 
the filter curves z = 0.5. The bins contain 7927, 14,750, 20,328, 
30,793, 19,648, and 8702 objects (from lower to higher redshifts). 
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Figure 16. Density plots of z vs. absolute R band magnitude of 
LRGs. Upper panel: all LRGs where Ha lies within the wavelength 
coverage of SDSS. Lower panel: Relative abundance of LRGs ful- 
filling logio(OIII/H/3) < 0.7 - 1.2(logio(NII/Ha) - 0.4) and hence 
showing star formation activity. 



emission hnes. Galaxies that fulfih log]^o(Oni/H/3) < 
0.7 - 1.2(logio(NII/Ha) - 0.4) are tagged as "star form- 
ing" . If they furthermore have an Ha equivalent width of 
equal or greater than 500 A, they are classified as "star- 
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For a definition of "star forming" see the text. Error bars indicate 
the Poissonian errors for the intermediate bins, and Binomial errors 
(with 84% confidence) for the outer bins which are not as highly- 
populated. The bins contain (from higher to lower luminosity) 1, 
6917, 23,429, 1751, 122, and 34 objects in total (star-forming and 
non-star-forming) . 
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Figure 18. Stacked rest frame spectra for S PSS LR Gs a t z < 0.1 
and within the same z — Mji bins as in Figures [131 [III andfTTl The 
SEDs were normalized to fg, interpolated to a resolution of 5 A , 
and weighted by the calculated fji (likewise as those in Figure |14| > 
in the stacking. Different colors represent different Mji bins. TTTe 
spectra are shown within the combined SDSS wavelength coverage 
for objects between z = 0.0 and 0.1, whereas plotting limits are 
the same as in Figure [14] SDSS filter transmissions are displayed 
by solid (z = 0.0) and dashed (z = 0.1) gray lines. 



burst". In our LRG sample only 10 objects fulfill the 
"starburst" criteria. In the following, we will refer to 
both tags ( "star forming" and "starburst") as star form- 
ing. Since the spectral wavelength coverage of SDSS is 
[3800 A, 9200 A], the Ha line drops out of the sample at 
z > 0.4. From the previously investigated ~ 140,000 
objects, ^ 120,000 remain in the sample when we cut 
off LRGs that have z > 0.4. From these, 4000 show 
indications of star formation. The upper panel of Fig- 
16 shows the density of LRGs in the z — Mr plane. 



and the lower panel displays the relative abundances of 
LRGs exhibiting star formation. Here, one can observe 
that the fraction of local star forming LRGs increases 
toward smaller luminosities. Additionally, we show in 



rors. The errors of the outer two luminosity bins, which 
do not contain as many objects, are Binomial with a 
confidence level of 84%. In principle, the highest and 
lowest Mr bins could be neglected, since they contain 
only 1 and 33 objects, respectively. Except for the re- 
gion of faint galaxies (the lowest luminosity bin), the 
frequencies of star-forming LRGs shows an increase to- 
ward fainter magnitudes. This is in agreement with what 
we would expect from Figure [T4j where we saw that the 
mean model SEDs for faint LRGs are bluer and show 
more star-forming activity, as well as with the results of 
the fi^ we performed in the last section. 
In addition to that we downloaded spectra for random 
subsamples of the objects with z < 0.1 and stacked them. 
Thereby we weighted them by the best fit SED flux in the 
r band as w e did with the superposition of our models in 
Section |7.2| In each luminosity bin the spectra of a sub- 
sample containing 100 objects (less, if the bin contains 
less objects) were superpositioned. We show the results 
in Figure [Tsj where we consider the added spectra in the 
rest frame wavelength range which is covered by SDSS 
for objects with z < 0.1. Comparing Figure [l8 to Fig- 



Figure 17 the percental fractions of star- forming LRGs 
for the same redshift-absolute R b and magnitude bins 

that we investigated in Section Figures [13] and [T4j The overall bias reaches {\Az\/{l + ^spec)) = 0.015. 
The errors in the intermediate bins are Poissonian er- We compared our results with photometric redshift 



ure[14, we observe that the stacked SDSS spectra show 
a very similar behavior as the average best fitting SEDs 
that were inspected previously, supporting that our novel 
SEDs are able to describe the data of this population. 
At greater redshift s {z > 0.1), no faint galaxies are ob- 
servable, and only 370 0.5 %) LRGs are classified as 
starforming. This region is therefore strongly dominated 
by objects that show no recent star formation. The vari- 
ations we saw in Figure 15 are not caused by different 
fractions of star forming galaxies (which are anyhow of 
very low abundance), but result from aging only, which 
was also support ed by the fit to BC03 models we per- 
formed in Sectior fr2l 

We also stacked spectra of objects from these z — Mr 
bins in the same manner as above. Unfortunately, one 
cannot see the effect of a flatter NIR and declining UV 
flux at low z compared to a declining NIR and increasing 
UV flux at high z as clearly as in Figure [T5| because the 
NIR part is only observable for objects with low and 
the UV only for objects with high z. However, the trends 
at low z (i.e., flat NIR flux) and at high z (increasing UV 
flux) are qualitatively there. 

8. SUMMARY 

In this work we analyzed the colors as a function of 
redshift of a number of publicly available red model 
SEDs in the SDSS filter system and compared them to 
the colors of spectroscopic SDSS LRGs at redshifts from 
z = to 0.5. We found that no investigated model set 
can describe the u — g^ g — f'i r — i, and i — z colors of the 
data at all redshifts. Thus we created new red model 
SEDs for five redshift bins on the basis of the LRG pho- 
tometry by fitting a superposition of a BC03 CSP model 
and a burst SED, while allowing for dust extinction of 
both components. We estimated photometric redshifts 
with the Bayesian template fitting code PhotoZ and a 
template set containing the created SEDs with according 
redshift and luminosity priors. The resulting catas- 
trophic outlier rate is 7^ = 0.12%, whereas the scatter 
reads aAz = 0.027 and ctnmad = 0.017 for non-outliers. 
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estimates available in the SDSS database, both for the 
template fitting approach and for the neural networks. 
Our (|Az|/(l -\- Zspec))^ (^Azi and (Tnmad are of the same 
order as the SDSS database template fitting values, 
although with a slightly higher catastrophic outlier rate. 
Nevertheless, the medium bias varies significantly less 
from = with our template SEDs. The total values 
of our scatter, bias and outlier rates marginally differ 
from the SDSS ANN(CC2) photometric redshifts and 
are worse than the results of ANN(Dl). However, the 
bias of ANN redshifts reaches higher values within some 
redshift ranges. 

We investigated the properties of the novel SEDs and 
demonstrated that they also describe the variety of 
GALEX NUV and FUV colors well. Performing an 
SED fit at the spectroscopic redshift of the LRGs, 
we observed that the flux weighted mean SEDs of 
local galaxies (z < 0.1) vary notably as a function of 
luminosity, meaning that faint galaxies show increasing 
indications of ongoing star formation. The stacked 
spectra from the SDSS show a very similar behavior. 
Furthermore, we find that the fraction of star forming 
SDSS LRGs at 2: < 0.1 increases with decreasing R band 
luminosity. This is in agreement with our finding that 
the mean SED is bluer for faint local galaxies. We also 
observe a difference in the flux weighted mean SEDs as 
a function of which is an effect of aging only and not 
due to recent star formation. In Appendices |A] and [B] 
we showed that BC03 models cover the data better than 
M09 and Mil in color space, which is why we favored 
the BC03 relative to the M09 and Mil model sets for 
the creation of our new SEDs. Finally, we also showed 
how our new template SED colors compare to the data 
in the ugri^ griz^ grri, and riiz planes. 
We would be happy to provide the photometric redshifts 
together with the rest frame magnitudes for the complete 
photometric LRG sample of SDSS on request. 
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APPENDIX 



sed fitting results with different model sets 

In Section [4] we created templates by fitting model SEDs to the photometry of the SDSS LRGs in order to produce 
SEDs matching the colors of the data. Here we want to show the fitting results and how well the models are able to 
cover the observed colors. In Figures |2] and [5] we saw that both the BC03 and the M09 model SEDs can represent 
the LRG data, at least for a large part of the investigated redshift range. Both sets have a large variety of ages (plus 
e-folding timescales and metallicities in the case of BC03 models), thus they offer the perspective to provide a well 
matching model for each redshift. 

In addition to that, iMaraston & Stromback ( 2Q11[ Mil) provide a large variety of single stellar populations (SSPs) 
created from different stellar libraries and with varying metallicities and ages. The Mil models are promising can- 
didates to represent the LRG data in terms of colors, wherefore we analyzed them in this regard and compared to 
the results with BC03 models. We could only use those Mil model SEDs that are provided for a wavelength range 
sufficiently large t o cover a ll SDSS filters at all LR G redshifts. These are two flux-calibrated empirical stellar spectral 

and MILES (Sanchez-Blazquez et al. 2006), and the theoretica lly derived MAR CS 

([2003| 



libraries, Pickles 



(Gustafsson et al. 



(| Pickle s||1998|) ; 
■11:^0081 tiEH?y. 



The Pickles models have solar metallicity, either an IMF of Chabrier 



or 



Salpeter, and were extended in the near-IR by Mil. A version with a theoretical UV part was additionally created 
by Mil. The MILES SSPs in question were extended in the near-IR and UV, have a revised IR slope, an d solar 
metalhcity. The MARCS models have metallicities of Z = 0.0004, 0.02, and 0.04, and IMFs of Salpeter and Kroupa| 
([2001) respectively... A more thorough description of the SSPs as well as the SEDs can be found in Mil and at 
pittp://www. icg.port .ac.uk/~inaraston/Mll/, These models are a subset of all Mil covering a wavelength range 
wide enough for photometric redshift estimation. 



We produced CSPs from the selected Mil models with the software EzGalF^ (Mancone fc Gon zalez||2012 ). We apply 
exponential SFRs and the same e-folding timescales that were used with BC 03 (Section 3T272). The whole range of 



of 

available ages, IMFs, and metallicities were exploited where possible in order to cover the greatest variety within the 
resulting SEDs. 

These CSPs were then fitted to the photometry in the same way as the BC03 CSPs in Section |4] (i.e., superposed by a 
burst model and with additional extinction) for each stellar library and IMF separately. The same was also done with 
M09 models. 

The SED fitting results for the different redshift bins and color spaces of M09/11 are shown together with those of 
the BC03 models in Figure 19 From the M09/11 models in all cases the one with the smallest result is chosen 
for the plot. The upper left panel of Figure \i9\ {z ^ 0.02) shows that the BC03 models cover the LRG colors better 
in this z range, with the exception of riiz^ wnere M09/M11 account for the lower right part that is not covered by 
BC03. Furthermore, the BC03 color contours trace the LRG distribution better than M09/11 in grri and griz. In 
the upper right and middle left panel the colors are shown for the redshift bins at z ~ 0.1, and z ~ 0.2. Both model 
sets perform similarly well, except for z ~ 0.2, where M09/11 produces a slightly wider spread in ugri. Apart from 
that, the BC03 contours follow well the density of the LRGs in the color spaces grri and griz. At z ~ 0.3 the M09/11 
models outperform BC03 in riiz^ whereas both describe the data colors similarly in the other spaces (middle right 
panel). Here, the M09/11 SEDs produce two peaks in riiz which cannot be observed for the LRGs. Finally, the lower 
panel of Figure [19] shows the redshift bin around 0.4, where M09/M11 still produce a wider spread in riiz^ but are not 
able to match the bluer part of the LRG'sr — i color. In addition to that, neither BC03 nor M09/11 models cover the 
bluer part of the i — z color. In Section 4.2 we already mentioned that the spread in the colors of the data cannot be 



completely accounted for by the photometric errors in this redshift range. The reason for the spread must thus be of 
physical origin. Unfortunately, not the whole range of i — z of LRGs at z ^ 0.4 can be fitted by the investigated models 
sets while simultaneously matching the g — r and r — i colors (where the latter has smaller errors than i — z). This is 
a feature of BC03 and M09/11 models, and the reason why we could not create SEDs with bluer i — z in Section [4] 
The distributions of the results of the fitting to the various model sets are shown in Figure 20 Each analyzed z 
bin is displayed in a different panel, whereas the distributions for all redshift intervals combined is shown in the 
lower right panel. BC03 outperforms the other models (M09 and Mil), regarding the values, in every redshift slice 
except for z ~ 0.4, where the Mil Pickles models fit the LRGs better, but only marginally. We would like to add that 
the Mil Pickles models fit the data best compared to the other M09/M11 SEDs, and outperform the Pickles models 
with a theoretical UV-part of the spectrum for each fitted object. 

Finally, the lower right panel of Figure 20 shows that the values are best also for all z bins combined, as expected 
from the previous results. Following these results, and because of the previous analysis of the colors of the produced 
SEDs, we adopted the BC03 models for the creation of photo-z template SEDs. 

COLOR COLOR RELATIONS OF THE NEW MODELS 



Here we present the color-color relations of our models and the LRG data. Figures [6j and 21 25 show the colors of 
the best fitting SED models, the preselected models, and the selected models after the removal of redundant SEDs in 
comparison to the LRG data itself in the ugri, griz, grri, and riiz color spaces. We display the LRGs with u band 
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Figure 19. Color distribution of LRGs (gray). Contours show the SED fitted BC03 (green), and M09/11 (red) models superposed with 
a burst model to fit the LRGs. 
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Figure 20. distribution of SEDfit results of different model-sets: BC03 (black), Mil MARCS stellar library with Salpeter IMF (blue), 
MARCS with Kroupa IMF (hght blue). Mil MILES (green). Mil Pickles (yellow). Pickles with theoretical UV (red), and M09 (light red). 
The different panels show the various redshift slices in which the fitting was performed. The lower right panel finally displays the results 
in all redshift intervals combined. 



Photo-zs of LRGs from the SDSS DR7 



19 




Figure 21. Color vs. color plots for SDSS data (blue) and their individual best fitting SEDf it-SEDs (red) for 0.0 < 2; < 0.04. Objects 
with u band errors lower than the median are indicated in light blue, whereas those with greater errors are dark blue. The mean error bars 
of all objects are shown in the lower right corners in dark blue, and the mean errors of objects with u band errors smaller than the median 
in light blue. The grid in r — i and the boundaries within which the models are selected are shown in black. The dark green dots are the 
preselected models, whereas the light green points represent the models that are left over after the removal of redundant SEDs. The yellow 
dots are models that shall account for objects outside the selected boundaries in ugri and griz. 

errors below the median and those with errors above with different colors to visualize the difference in the populated 
loci within the ugri space, which becomes evident at z ^ 0.3 (Figure 24, upper panel) and z ^ 0.4 (Figure [25] upper 
panel). 

The grids we imposed on the ugri plane within which we selected templates (see Section 4.2) are displayed in the 
upper left panels of Figures [6], and 22 -25, whereas the color boundaries in ^ — r and i — z are drawn in the lower panel 



of Figure [6| in Figure [21] and in theTower left and in the right panels of Figures [22 -25 
In order to account formDjects with peculiar colors which lie outside the boundaries as well, we select additional SEDs 
by eye and confirm that they yield a photometric redshift bias lower than 0.01, when used as single template for the 
subcatalog from the same z range. If they do not meet that requirement, we chose the SED of another object from 
the color neighborhood of the first one, test for the bias, and so on until we find a model from that color locus which 
creates only a small Zphot-bias 



Bright cluster galaxies (BCGs) with very low UV-excess (i.e., in the lower part of the UV red sequence; see Section 
Figure 12) at z ~ 0.1 are not matched by the models already selected for that z range. But we find that two S^ 
originatmg from z ^ 0.2 are able to cover the considered region. Therefore, we add these templates to z ^ 0.1, 
we assign the same redshift priors as for z ~ 0.1 models in addition to their former priors (see Section [2|, thereby 
covering the lower UV red sequence. The colors of these two models are shown by yellow triangles in Figure [22] 
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Figure 22. Color vs. color plots for SDSS data (blue) and their individual best fitting SEDf it-SEDs (red) for 0.08 < z < 0.12. Objects 
with u band errors lower than the median are indicated in light blue, whereas those with greater errors are dark blue. The mean error bars 
of all objects are shown in the lower right corners in dark blue, and the mean errors of objects with u band errors smaller than the median 
in light blue. The grid (upper panel) and the boundary (lower panel) within which the models are selected are shown in black. The dark 
green dots are the preselected models, whereas the light green points represent the models that are left over after the removal of redundant 
SEDs. The yellow dots are models that shall account for objects outside the selected boundaries in u qri and griz. Yellow triangles are 
models from z ~ 0.02 that lie in the lower part of the UV red sequence of R12 (see Section [TT] Figure|T2j and are added additionally in 
this redshift range. 




Figure 23. Color vs. color plots for SDSS data (blue) and their individual best fitting SEDf it-SEDs (red) for 0.18 < z < 0.22. Objects 
with u band errors lower than the median are indicated in light blue, whereas those with greater errors are dark blue. The mean error bars 
of all objects are shown in the lower right corners in dark blue, and the mean errors of objects with u band errors smaller than the median 
in light blue. The grid (upper panel) and the boundary (lower panel) within which the models are selected are shown in black. The dark 
green dots are the preselected models, whereas the light green points represent the models that are left over after the removal of redundant 
SEDs. 
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Figure 24. Color vs. color plots for SDSS data (blue) and their individual best fitting SEDf it-SEDs (red) for 0.28 < z < 0.32. Objects 
with u band errors lower than the median are indicated in light blue, whereas those with greater errors are dark blue. The mean error bars 
of all objects are shown in the lower right corners in dark blue, and the mean errors of objects with u band errors smaller than the median 
in light blue. The grid (upper panel) and the boundary (lower panel) within which the models are selected are shown in black. The dark 
green dots are the preselected models, whereas the light green points represent the models that are left over after the removal of redundant 
SEDs. 




Figure 25. Color vs. color plots for SDSS data (blue) and their individual best fitting SEDf it-SEDs (red) for 0.38 < z < 0.42. Objects 
with u band errors lower than the median are indicated in light blue, whereas those with greater errors are dark blue. The mean error bars 
of all objects are shown in the lower right corners in dark blue, and the mean errors of objects with u band errors smaller than the median 
in light blue. The grid (upper panel) and the boundary (lower panel) within which the models are selected are shown in black. The dark 
green dots are the preselected models, whereas the light green points represent the models that are left over after the removal of redundant 
SEDs. 



